Quantum metrological limits via a variational approach 
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The minimum achievable statistical error in the estimation of physical parameters is determined 
by the quantum Fisher information. Its computation for noisy systems is still a challenging problem. 
Using a variational approach, we present an equation for obtaining the quantum Fisher information, 
which has an explicit dependence on the mathematical description of the noise. This method is 
applied to obtain a useful analytical bound to the quantum precision in the estimation of phase- 
shifts under phase diffusion, which shows that the estimation error cannot be smaller than a noise- 
dependent constant. 
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Quantum metrology QJH] deals with the estimation of 
parameters taking into account the constraints imposed 
by quantum laws. The estimation is based on measure- 
ments made on probe systems undergoing a parameter- 
dependent process. For a given measurement scheme, 
the uncertainty in the estimation of a parameter is lim- 
ited by the Cramer-Rao bound, which is proportional to 
the inverse of the square root of the so-called Fisher in- 
formation (FI) [IHZ]. The maximization of FI over all 
measurement strategies allowed by quantum mechanics 
leads to a non-trivial quantity: the quantum Fisher in- 
formation (QFI). The determination of this quantity is 
central to quantum metrology. It allows, for instance, 
the establishment of ideal benchmarks for the statistical 
error in the estimation of parameters, which can be used 
by experimentalists to evaluate the performance of a real 
experiment. A systematic approach to calculate the QFI, 
using the symmetric logarithmic derivative (SLD) oper- 
ator [H [2] , was developed in Ref . [3] . This approach has 
allowed large advances on quantum metrology [5] IS] • For 
unitary processes, it leads to simple analytical expres- 
sions. This is not the case, however, for noisy processes, 
which often require numerical calculations. 

When the unknown parameter is associated with a 
physical process, the ultimate limit for the precision in 
its estimation is given by a further maximization of the 
QFI over all initial probe states (given some constraint, 
e. g. a fixed average energy). These two maximizations 
make the determination of that ultimate limit a labo- 
rious numerical task. Recently, an alternative to solve 
this problem was presented in Ref. [lOj : given a mathe- 
matical description of the quantum parameter-dependent 
process by a set of Kraus operators [IT] , an upper bound 
to the QFI can be calculated; the true value of the QFI is 
obtained by minimizing this upper bound over all equiv- 
alent Kraus representations of the process. 

In this article, we present a variational approach, based 
on purification techniques [T2], to calculate the QFI 
through the minimization of upper bounds. These upper 
bounds correspond to the QFI associated with pure states 
in the enlarged Hilbert-space of all purifications of the 



original probe state. An important advantage of our ap- 
proach is that it results in a general prescription for per- 
forming that minimization: beginning with an arbitrary 
purification of the probe state, the optimum purification 
that minimizes these upper bounds and leads to the QFI 
can be found through the solution of a Sylvester equation 
|13) . Such a prescription, which leads to an alternative 
way of expressing the quantum Fisher information, actu- 
ally solves the minimization problem posed in Ref. [lOj . 
As a concrete example, we use this approach to delimit 
the ultimate precision bounds on the estimation of phase- 
shifts in the presence of phase-diffusion. This problem 
was addressed recently in [14], where these bounds were 
found numerically for initial Gaussian probe states, and 
posteriorly confirmed experimentally in the special case 
of initial coherent states [TS]. Here, we show an analyt- 
ical and non-trivial lower bound for this quantum limit, 
which is valid for any probe state. This bound reveals 
a drastic effect of phase-diffusion on phase estimation: 
the accuracy, even though dependent on the energy of 
the probe state, cannot be better than a noise-dependent 
constant. 

The estimation of a parameter is based on experimen- 
tal results of measurements on a quantum probe state, 
after it has been submitted to a physical operation that 
depends on the value of the parameter. Usually, it is 
previously assumed that the possible values of the pa- 
rameter to be estimated, denoted here by x, lie within 
a certain continuous interval. It is also assumed that 
one knows the precise dependence of the physical op- 
eration on the value of the parameter x, so that, for 
a given initial probe state, one knows the transformed 
state p(x). Finally, one assumes a specific measurement 
device, mathematically represented by positive operator- 
valued measures (POVM) {Eh} [12]. The goal, therefore, 
is to obtain an accurate estimation, x es t(k), of the true 
value Xtrue of the unknown parameter, from a set of ex- 
perimental results k, using a given rule. The variance 
between an estimate and any possible value of the pa- 
rameter is Sx = yj ((x cst — x) 2 ), where (•) = ^2 k Pk(x)* 
and Pk(x) = Tr[p(x)Ef.]. Here, pk(x) is the probabil- 
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ity of obtaining the experimental result k given that 
the parameter value is x. This variance may be con- 
sidered a merit quantifier for the error estimation as 
a function of x. For any unbiased estimation (i. e. 
(x cs t) = x), the statistical error is limited by the Cramer- 
Rao bound [5] [7]. For an experiment with v repeti- 
tions, this bound is given by 6x > l/y/ vF(x), where 
F( x ) = J2kPk( x )[d \n[Pk(x)]/dx] 2 is the FI. Under very 
general assumptions [SHI] , it can be shown that it is pos- 
sible to saturate the Cramer-Rao bound, at least in the 
asymptotic regime (v —> oo). 

The QFI is defined by the maximum of the FI over 
all possible measurement strategies allowed by quantum 
physics: 

T Q [p(x)]=nmxF[p(x);{E k }] . 

{E k } 

The respective quantum version of the Cramer-Rao in- 
equality [IB], 5x > 1/ \JvFq(x), settles, therefore, a limit 
to the statistical error that cannot be overcome by any 
strategy of estimation, for given physical process and 
quantum probe. 

If the transformed probe state is pure, p(x) = 
\tl)(x)){ij}(x)\, the correspondent expression of the QFI 
is[THirjI]: 
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If the state p(x) is not pure, the SLD approach does 
not lead in general to such a simple analytical expres- 
sion. On the other hand, it is always possible to en- 
large the size of the original Hilbert space S and build 
a pure state \$s,e(x)) m the enlarged space S + E 
that fulfills the condition TtePs.e(x) — ps{x), where 
PS,e(x) — \®s,e(x))(§s,e(x)\, and the trace is taken only 
on the S-Hilbert space [12] • We have added the label S 
to the state p{x) of the system, in order to distinguish it 
from states in space S + E. The state \&s,e(x)) is called 
a purification of ps(x). The space E may be interpreted 
as an environment for system S. 

Because taking the trace over E may be interpreted as 
discarding information on part of the total space S+E, a 
physically motivated upper bound Cq of J~q[ps(x)] can 
be obtained: 

Cq [Ps,e{x)\ = FqIPsM 3 *)] ^ Fqlpsix)] , 

this inequality being valid for any purification of ps{x). 
Physically, this is due to the fact that when a system plus 
an environment are monitored together, the information 
obtained about an unknown parameter cannot be smaller 
than the information acquired when only the system is 
measured. Since Cq depends on the purification chosen, 
the best upper bound that can be obtained with this 
strategy is given by the minimum of Cq over all possible 



purifications of p(x). In the supplementary material of 
Ref. [10] . it is shown that this minimum is equal to the 
QFI: 



{\$s,e(x))} 



if the dimension of E is at least equal to the dimension 
of S. 

It is possible to determine the value of QFI by minimiz- 
ing the upper bound Cq over all purifications of ps(x) 
[TU] , This is, in general, a challenging task and a con- 
crete prescription to do it would be welcome |19) . In the 
following we present such a prescription. Our procedure 
starts by establishing a relation between all purifications. 
As shown in Ref. [12], there is always a unitary oper- 
ator ue{x) in i?-space that connects two purifications 
\^S,e(x)) and \&s,e(x)) of the same state ps{x): 

\$s,e{x)) = ue(x)\$s,e(x)) ■ 

Therefore, given a purification \<&s,e(x)), the QFI may 
be found by minimizing Cq[ue(x) ps,e(x)u e (x)} over all 
unitary operators on _E-space. The physical role of ue(x) 
is to erase all non-redundant information about the pa- 
rameter x that has been leaked from the system into the 
larger space S + E. 

It will be useful to define a Hermitian operator h(x), 
which acts only in the i?-space, by 

h{x)=i y ' ue(x) , 

and another Hermitian operator H(x), which acts in the 
whole S + S-space, by 

.d|*s,E(a:)> 
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Using the definitions above, we may write Cq as 
Cq=4([H(x)~(H(x))z} 2 )<, , 
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where H(x) = H(x) — h(x), and the averages are taken 
over |$ s,e{ x )) ■ Eq. (|2j) expresses the quantum Fisher in- 
formation is terms of the variance of the generator of the 
parameter-dependent unitary evolution in the enlarged 
space. From Eq. ([2]), we conclude that the minimiza- 
tion of Cq over all unitary operators ue(x) is equiva- 
lent to the minimization of Cq over all Hermitian oper- 
ators h(x) that act on _E-space. This minimization is 
a mathematical optimization problem in semi definite 
positive quadratic programing, which can be efficiently 
solved numerically [20], since the operator h{x) appears 
as a quadratic function in Cq. Thereupon, it is possible 
to find an equation for the optimum Hermitian opera- 
tor h pt{x) that minimizes Cq. Taking, without loss of 
generality, (h op t(x))<£ — (H(x)}$, one finds that: 



h pt(x)p E (x) + PE{x)h opt (x) 



Tr s {V[ps, E (x)]} , 



(3) 
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where pe{%) = Trs[|<i>s,E(x))(<i>s,E(x)|] is the reduced 
density matrix in the I?-space, and T> [ps,e(x)] is defined 
as: 



Markov limit, by 
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Equation ([3]), when expressed in terms of the matrices 
associated to the corresponding operators, is a Sylvester 
equation [13j . It depends only on the degrees of freedom 
of E. After determining h opt (x), the QFI may be finally 
expressed as 



T Q [ps{x)] = C q [ps,e(x)} - 4(A^ pt (*)> 



This is a novel expression for the quantum Fisher infor- 
mation, which relates it directly to the QFI correspond- 
ing to an unitary evolution of the enlarged system, and 
shows that the non-redundant information in \<&s,e{x)) 
about the parameter x is given by four times the vari- 
ance of h op t(x). As compared to the expression for the 
QFI of the system alone, given by the SLD approach, it 
displays explicitcly the mathematical description of the 
noise process, through the purification \&s.e{x)}- 

The usefulness of this method is not restricted to ob- 
taining an equation for the exact evaluation of QFI. Since 
our approach is variational, whenever it is too hard to 
find a solution of Eq. ^ , we may still obtain satisfactory 
and non-trivial analytical upper bounds to QFI. Indeed, 
based on (Ji), one may guess approximations for h opt (x) 
that depend on variational parameters, so that the min- 
imization is made on sub-classes of operators h{x). This 
procedure leads to bounds for QFI, given by the min- 
ima of § over these sub-classes. It also allows an iter- 
ation procedure for getting progressively better approxi- 
mations to the QFI of the system. 

The estimation of phase-shifts is a central problem 
in quantum optics, metrology and quantum communi- 
cation [T4l IT51 [2"T] . It is important, for example, in the 
use of light interferometers as part of detectors of gravita- 
tional waves [22] • When planning such experiments, it is 
essential to take into account the unavoidable influence 
of noise on the ultimate precision limits. Besides pho- 
ton losses, phase diffusion is another relevant source of 
noise in optical phase measurements and must be taken 
into account. For incoming Gaussian states, a numerical 
study of the effect of phase diffusion on the ultimate limit 
of precision for phase estimation was presented in [13] . 

In the following we apply our approach to this problem 
in order to derive non-trivial analytical lower bounds to 
this limit, which are valid for any probe state. For con- 
creteness, consider an initially pure probe state ps — 
\i ) s){'4 } s\, corresponding to a generic harmonic oscillator, 
which undergoes a phase-shift <f> due to some physical 
process. The resulting state of the probe, in the pres- 
ence of phase diffusion noise, may be represented, in the 
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where p m ,n is the matrix element of the initial probe 
state in the Fock basis of the system and j3 quantifies the 
degree of diffusion present in the process (from /3 = 0, 
corresponding to no diffusion, to /3 = oo corresponding 
to maximum diffusion). In order to obtain a possible 
purification of ps(4>) based on physical insight, we con- 
sider an interferometer with a dispersive plate producing 
a difference of phase 4> between its two arms, and model 
the phase-diffusion of the initial probe state through the 
effect of the radiation pressure on one of the interferom- 
eter mirrors. The interaction between the light field and 
the mirror is taken as proportional to uxe, where h is 
the photon number operator and xe is the dimensionless 
position operator of the mirror. In this model, the final 
state of the combined system of the probe and the mirror 
is given by 
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where \0e) is the initial state of the mirror, which we 
assume to be the ground state of a quantum oscillator. 
It is straightforward to see that this state is, indeed, a 
purification of ps{4>)- The value of Cq[ps,e(4 > )] mav be 
now calculated directly through Eq. (nl: Cq — 4An 2 , 
where An 2 is the variance of the photon number operator 
in the initial probe state. Notice that this purification 
leads to a trivial upper bound to the QFI of p(4>) as 
it is equal to the QFI in the absence of phase diffusion 
(0 = 0). 

In order to obtain a tighter upper bound to the QFI 
of p((f>), we consider possible approximations to h opt (x) 
by analyzing more closely Eq. ([3]). The reduced den- 
sity matrix of the mirror associated with the purification 
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where \i\2f3n) E is a coherent state with amplitude 
\/2/3n. The right-hand side of Eq. (B is Tr s D[p S|B ] = 



[ib E /(2V2/3)}p E - ip E [b E /(2V2(3)} 7 where b E = (x E + 
iPe)/V2, with being the dimensionless momentum 
operator of the mirror. Notice that the solution of Eq. ^ 
for h opt (x) would be trivial if ibE were an hermitian op- 
erator. However, in the asymptotic regime, y/2(3n ^> 1, 
the operator Pe/(2/3) when applied to pe produces a re- 
sult quite similar to that of ib/(/3y/2). So, we may guess 
that u E (<t>; A) = e -^W(2,3) ( with ^ being a var i a ti nal- 
parameter) would be a reasonable candidate to erase part 
of non-redundant information in \^s,e(<I>))- in this case 
the upper bound of QFI is: 
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FIG. 1: Comparison between upper bound Cg ax and the max- 
imum quantum Fisher information J-q^ in [14] as a function 
of the average number of photons N . The dots stand for 
the values obtained in [T3], while the full lines correspond to 
C™. The inset displays the two quantities up to TV = 30, 
which was the range considered in [14] . From top to bottom, 
P 2 = 0; 5 x 10~ 6 ; 5 x 10~ 5 ; 5 x 10 -4 . 



The optimal value of A that minimizes Cq is A opt = 
8An 2 /3 2 /(l + 8An 2 /? 2 ). Then, taking the inverse of the 
square root of Cg pt , one gets a non-trivial bound for 
the precision of phase-shift estimation in the presence 
of phase-diffusion, valid for any input state: 

H pd >^l/{vC^) = ^l . (4) 

This inequality shows that the error in this estimation is 
limited by a well-known formula for independent errors, 
displaying clearly the effects of the intrinsic probabilis- 
tic feature of quantum mechanics, l/(4An 2 ), and of the 
genuine phase-diffusion noise, 2/3 2 . 

An important property of the bound shown above is 
the presence of a constant term. This means that the 
presence of phase diffusion is, in general, more detrimen- 
tal to phase-shift estimation than the presence of photon- 
losses, when the error goes to zero as the average number 
of photons goes to infinity. 

From Q, it follows that, for Gaussian states, one may 
obtain a bound that depends explicitly on the average 
photon number N: 



2/3 2 + 



87V(iV + 1) 



since, for these states, An 2 < 2N(N + 1). 

We compare in Fig. [T] Cg lax with the maximum quan- 
tum Fisher information J r g ax obtained numerically in 
[M] for the best Gaussian probe states with given av- 
erage photon number N. Good qualitative and quantita- 
tive agreements between them is observed for N up to 30, 
which is the range considered in [TJ . On the other hand, 
our bound, being analytical, allows one to obtain a better 



insight, for any value of iV, on the ultimate limit for phase 
estimation in the presence of phase diffusion. In partic- 
ular, the saturation of our bound when 16(iV/3) 2 ^> 1 is 
clearly displayed. 

In conclusion, we have presented in this article a varia- 
tional method to determine the quantum Fisher informa- 
tion by minimizing upper bounds to this quantity, and 
have given a general prescription to perform this mini- 
mization. It bears all the advantages of variational meth- 
ods, which lead to useful analytical bounds in situations 
where an exact solution cannot be found analytically. 
We have applied this method to phase-shift estimation in 
the presence of phase-diffusion and have obtained a non- 
trivial lower bound to its statistical error. This bound, 
which agrees with and goes beyond published numeri- 
cal results, shows that there exists a constant limit to 
this error, which depends only on the strength of the 
phase diffusion. We believe that the method proposed 
here might be very useful to determine the fundamental 
precision limits in quantum metrology in the presence of 
noise. 
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